Azimuthal spin asymmetries in light-cone constituent quark models 
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I. INTRODUCTION 

The composite nature of the nucleon has been explored for a long time by means of deep inelastic scattering (DIS) 
_ of a lepton beam in the Bjorken regime, i.e., when q and P denote the four-momentum transfer and the nucleon 
II . momentum, in the limit of P ■ q and = —q^ oo while x = Q^/ {2P ■ q) is fixed. As a consequence of the high scale 
CIh' Q, scattering occurs in a collinear configuration between the incident lepton and a single 'parton' in the nucleon. The 
1^ , factorization theorem allows the inclusive DIS cross section to be expressed as a convolution of two contributions: 
one corresponds to the hard process occurring at short distance between probe and parton; the other accounts for 
the coherent long-distance interactions between parton and target, and is described in terms of parton distributions. 
T— I ■ At leading order (leading twist) x can be interpreted as the fraction of the longitudinal momentum of the parent 
^ ' (fast-moving) nucleon carried by the active parton, and one may distinguish three kinds of parton distributions. Two 
I of them are well-known from measurements of structure functions in DIS and other processes: fi{x) is the number 
' density of unpolarized partons with longitudinal momentum fraction x in an unpolarized nucleon, and gi{x) gives 
, the net helicity of partons in a longitudinally polarized nucleon. The third one, the (chiral-odd) transversity hi(x) 
• ' describing the number density of partons with polarization parallel to that of a transversely polarized nucleon minus 
,; the number density of partons with antiparallel polarization, requires a quark helicity flip that cannot be achieved in 
• the inclusive DIS. Other processes have to be explored for that. 

I However, in addition to the information on the longitudinal behaviour in momentum space along the direction 
• • ■ in which the nucleon is moving, a complete three-dimensional picture of the nucleon also requires knowledge of 
the transverse motion of partons [1, 2]. A full account of the orbital motion, which is also an important issue 
to understand the spin structure of the nucleon, can be given in terms of transverse-momentum dependent parton 
^ ' distribution functions (TMDs). There are eight leading-twist TMDs fi{x,pT), f^{x,pT), giL{x,pT), 9it{x,pt), 
9^ 1 hi{x,pT), hj-j^{x,pT), hj^rpixjpx), hj-{x,pT) [3]. Two of them, the Boer-Mulders and Sivers functions hi{x,pT) and 
fij'{x,pT) [3, 4], are T-odd, i.e. they change sign under nai've time reversal, which is defined as usual time reversal, 
but without interchange of initial and final states. The other six leading-twist TMDs are T-evcn. 

In order to be sensitive to intrinsic transverse parton momenta it is necessary to measure adequate transverse 
momenta of the produced hadrons in the final state, e.g., in processes like semi-inclusive lepton-nucleon DIS (SIDIS), 
hadron production in e^e^ annihilation or the Drell-Yan processes in hadron-hadron collisions [1-21]. 

Here, factorization has been proved at leading twist [22-24] allowing to access information on TMDs as well as 
on fragmentation functions (FFs) describing the hadronization process of the hit quark decaying into the detected 
hadrons. At leading twist, the fragmentation of unpolarized hadrons is described in terms of two fragmentation 
functions, Di{z, Kt) and H^{z, Kt), where z is the energy fraction taken out by the detected hadron and Kt ~ \Kt\ 
its transverse momentum. The function Di{z, Kt) describes the decay of an unpolarized quark, whereas the Collins 
function H^{z, Kt) describes a left-right asymmetry in the decay of a transversely polarized quark [7-9]. 

By measuring the angular distribution of produced hadrons, in SIDIS it is possible to access information on all eight 
leading-twist TMDs in combinations with the two leading-twist FFs. Restricting ourselves to the one-photon-exchange 
approximation and considering spin degrees of freedom such as the beam helicity and the target spin, the contraction 
between the lepton and hadron tensors in the SIDIS lepton-nucleon cross section can be decomposed in a model- 
independent way in terms of eighteen structure functions, thus exhibiting a non trivial azimuthal dependence of the 
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detected hadron around the (space-like) direction defined by the virtual photon [10, 25-29]. According to factorization 
each of the leading-twist structure functions can be conceived as a convolution between one TMD and one FF. Since 
structure functions enter the cross section with a defined angular coefficient, they can be accessed by looking at specific 
azimuthal SIDIS asymmetries. This has become now a powerful tool for studying the three-dimensional structure of 
the nucleon [30-47], and many more data are expected to come in the future. The remarkable experimental progress 
was accompanied by and motivated numerous theoretical and phenomenological studies in literature [48-84] . 

In this paper we compute azimuthal spin asymmetries due to the T-even transverse-momentum dependent parton 
distributions functions. For that we use the predictions from the light-cone constituent quark model (CQM) of 
Rcf. [77], which has been successfully applied also in the calculation of electroweak properties of the nucleon [85] 
and generalized parton distributions [86]. Such a model, based on the light-cone wave function (LCWF) overlap 
representation of TMDs, is well suited to illustrate the relevance of the different orbital angular momentum components 
of the nucleon wave function for the respective observables. To best of our knowledge, this is the first attempt to 
describe all leading-twist spin asymmetries in SIDIS due to T-even TMDs in a single approach. We include also 
studies of the far better known coUinear double spin asymmetries Ai and All — not only for sake of completeness, 
but also to demonstrate the capability of the approach to describe reliably the gross features of spin effects in the 
nucleon. 

The article is organized as follows. In Sec. II the relevant definitions of azimuthal asymmetries in SIDIS are recalled. 
In Sec. Ill the main ingredients of the light-cone CQM of Ref. [77] are reviewed, and the results for T-even TMDs 
discussed. In Sec. IV we first study the collinear double spin asymmetries Ai and All- This and the following Sec. V 
devoted to a discussion of the transverse-momentum dependence of the TMDs, help to assert the range of applicability 
of the approach. In Sees. VI-IX we evaluate the leading-twist azimuthal double and single spin asymmetries due to 
T-even TMDs — focusing on their x-dependence. Sec. X exemplifies how the approach can be applied to make 
predictions for the transverse hadron momentum dependence of spin asymmetries. Concluding remarks are given 
in Sec. XI. Finally, a more detailed discussion about the model predictions for the sensitivity of the azimuthal 
asymmetries on different orbital angular momentum components are discussed in the Appendix. 



II. SPIN AND AZIMUTHAL ASYMMETRIES IN SIDIS 



The SIDIS process is sketched in Fig. 1. Let us denote 
the momenta of the target, incoming and outgoing lepton 
by P, I and I' and introduce the four-momentum transfer 
q ^ I - I' with = Then the relevant SIDIS vari- 

ables are defined as a; = Q'^/{2P ■ q), y ~ {P ■ q)/ {P ■ I) and 
z = {P ■ Ph)/{P ■ q)- The component of the momentum of 
the produced hadron transverse with respect to the virtual 
photon is denoted by Phi. and Ph± = 1-P?i±l- 

The SIDIS cross section (differential in x, y, z and the 
azimuthal angle (\)h of the produced hadron defined in Fig. 1) 
has the following general decomposition [10, 26], where ctq 
is the spin- and ^-independent part of the cross section, and 
where the dots indicate power suppressed ('subleading-twist') 
terms. 



HADRON PRODUCTION PLANE 




LEPTON SCA TTERING PLANE 



FIG. 1: Kinematics of the SIDIS process IN I'hX and 
the definitions of azimuthal angles in the lab frame. 
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In A'^Y^^ the index X describes the beam polarization, which is unpolarized (U) or longitudinal (L, characterized 
then by the beam helicity A). The index Y denotes the target polarization, which is unpolarized (U), longitudinal (L) 
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or transverse (T) with respect to the virtual photon. In experiments the target is polarized with respect to the beam, 
of course, but this is up to corrections of 0{1/Q) the same. As we shall deal with leading-twist observables, such 
corrections will be neglected through out. The superscript 'weight' reminds us of the kind of angular distribution 
of the produced hadrons with no index indicating an isotropic (/)-distribution, and 05 is the azimuthal angle of the 
target's transverse polarization vector, see Fig. 1. The asymmetries are defined in terms of structure functions, 
^weight _ ^wcight^^^^^ latter have the following partonic (tree-level) description in the Bjorkcn-limit [3, 11] 
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where h± = Pfi±/Ph± and M (m/i) is the mass of the nucleon (produced hadron). The convolution is defined as 



w j J 



= J d^ppj d^Kr S^^\zpp + Kp-P^^)w{pp, KT)Y,elxf{x,pT) ^iz^Kr) 



(11) 



where pr = \Pt\- These convolution integrals can be solved analytically only in the case of the structure functions 
Fuu and Fll, in which case the weight function w in Eq. (11) is simply unity. The unpolarized cross section is given 
in terms of Fuij by 



dxdydzd(l)h \^ ^ ' 2 
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Fuuix,z) , Fuu{x,z)=J2elxfnx)D'i{z) 



Recalling the notation gi{x) = J d^pp gip{x,pT), the double spin asymmetry All is given as follows 

Fll _ EaelxgmDii^) 



All = 



Fuu Y.a^lxmx)Dl{z) 



(12) 



(13) 



If no hadron is observed in the final state, the inclusive version of the double spin asymmetry (13) is commonly 
referred to as Ai and given by 



Eg elxgl{2 



(14) 



For all other structure functions (5-10) the convolution integrals cannot be solved analytically, unless one assumes 
models for the transverse parton momentum dependence of TMDs. A popular model is the Gaussian Ansatz, where 
one assumes 



fix,PT)=rix) ' r{z,KT) = r{z) 

1T{p^p{])) TT {K^{J)) 



(15) 



for some generic transverse parton momentum dependent distribution j'^{x,prp) and fragmentation J°'{z,Krp) func- 
tions. This is, of course, a crude approximation. However, besides being convenient [11], this Ansatz is also phe- 
nomenologically useful, provided the transverse hadron momenta are small compared to the relevant hard scale. 
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{Ph±) <C Q in SIDIS, and one is interested in catching the gross features of the effects [52]. A high precision de- 
scription of px-cSects requires methods along the QCD-based formahsm of [6], see [87] and references therein for 
examples. 

Using this Ansatz we obtain the following results 

^.i„(0.+0s) ^ B,Y.elxK{x)Ht^"'^'^{z), = (18) 
^UU - ^2 2_^e^xh^ {x)H^ (zj, B2- 1 + R{h^) ' ^ ' 



psin(30fe-0s) _ p, 2 ^ . -L(l)a/ \ T:r-L(l/2)a/ s p _ 3 

^^t/T --Bal^e.xV (-)^i ^^-{i?(/,i-^)i/3 + i?(/,i-^)-i/3}3/2' (21) 



where 

f{x)= j i'vrfii.PT). i'^"{i) = jifpr^fi^.-KT). (22) 

''."'"°(-'' - / '-'^^ '^") ^ ^ 

Of course, other models of pr-dependence can also be assumed. In that case, however, one typically cannot solve the 
convolution analytically, as in Eqs. (16)~(21), and has to use numerical integration. 

Notice that one could avoid the model dependence at this point by including adequate powers of transverse hadron 
momentum in the weights of the asymmetries [3]. The analysis of such P;i_L-weighted asymmetries is more involved, 
and so far only preliminary data not corrected for acceptance effects have been shown [38]. 



III. TMDS IN THE LIGHT-CONE CONSTITUENT QUARK MODEL 

A convenient way to describe parton distributions is to use the representation in terms of overlaps of LCWFs. This 
representation can be viewed as a generalization of the famous Drell-Yan formula for electromagnetic form factors [88] , 
and it was recently derived and applied in phenomenological calculations for the generalized parton distributions [89, 
90] and transverse-momentum dependent parton distributions [77, 91]. In practice, this representation becomes useful 
in phenomenological applications where one can reasonably truncate the expansion of the hadron state to the Fock 
components with a few partons. In our approach, we consider the minimum Fock sector with just three valence 
quarks. This truncation allows to describe the parton distributions in those kinematical regions where the valence 
degrees of freedom are effective, while the contributions from quarks and gluons are suppressed. 

The three-quark component of the nucleon has been studied extensively in the literature [92-98] in terms of quark 
distribution amplitudes defined as hadron-to-vacuum transition matrix elements of non-local gauge-invariant light- 
cone operators. Unlike these works, the authors of Refs. [91, 99, 100] considered the wave-function amplitudes 
keeping full transverse-momentum dependence of partons and proposed a systematic way to enumerate independent 
amplitudes of a LCWF given a particular parton combination. Within this general classification scheme, one finds 
that the nucleon state with three valence quarks has six independent scalar amplitudes which serve to parametrize 
the contribution from the four different orbital angular momentum components compatible with total angular 
momentum conservation, i.e. = 0, ±1,2. An application of this method has been developed in Ref. [77] for the 
calculation of the TMDs within a light-cone CQM which will be used here to make quantitative estimates of the 
azimuthal asymmetries. The key ingredient of the model is to derive the LCWF by boosting equal-time model wave 
function. The equal-time wave function is constructed as a product of a momentum wave function which is in a pure 
S-wave state and invariant under permutations, and a spin-isospin wave function which is uniquely determined by 
SU(6) symmetry requirements. The corresponding solution in light-cone dynamics is obtained through the unitary 
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Mclosh rotations acting on the spin of the individual quarks. By applying the Melosh rotations, the Pauli spinors of 
the quarks in the nucleon rest frame are converted into light-cone spinors. The effects of the relativistic spin-dynamics 
are evident in the presence of spin-flip terms in the Melosh rotations generating non-zero orbital angular momentum 
components which can be mapped out into six independent scalar amplitudes. The explicit expressions of these 
light-cone amplitudes can be found in Ref. [77] , while the corresponding results for the TMDs are given by 

/f(x,pT) = N^' f d[X] S{x-x,)S{pt-p^,) mx,},{p^,})\', (24) 



(25) 



gtrix^PT) = P'' J d[X]Six-X3)5{pT-p^3) 



2M{m + xMo) 2 



(26) 



h1{x,PT) = P'^ [ d[X] S{x-x,)5{pr-p^s) I ^r^AJ^°r 2 m{^.},{P±M'. (27) 

/O A/f2 
d[X] 5{x-x,)6{pr-p^,) ^^^^^m.Y^pI (28) 

hi£{x,PT) = -P^ I d[X] 5{x-x,)5{pr-p^,) '^^''i^Zi^2^1'l m{^^}AP± .})\', (29) 
where we introduced the integration measure 

dm = dx,dx.d.35 (^i-E-^) ^^^^^^^^f^s (tp^^ ■ (30) 

In Eqs. (24)-(29), Mq is the mass of the non-interacting three-quark system, and m the constituent quark mass. 
Furthermore, the flavor dependence is given by the factors = 2, N'^ = 1, and = |, P'' = — i, as dictated 
by SU(6) symmetry. A further consequence of the assumed SU(6) symmetry is the factorization in Eqs. (24)- (29) 
of the momentum-dependent wave function 'ip{{xi}, {p^.^}) from the spin-dependent factor arising from the Melosh 
rotations. Thanks to this factorized form one finds the following relations 

2/j?(x,pt) = 5?L(a^,PT) + ^/r(x,pT), (31) 

^/f(x,pT) = K{x,pT)-^hi^-{x,pTl (32) 
hilix.pr) = ~giT{x,PT)- (33) 

Eq. (31) is a generalization of analogous relations discussed in [55, 101] and was also rcderivcd together with Eq. (32) 
in Ref. [76]. Eq. (33) was already found in the diquark spectator model of Ref. [75]. In QCD the various TMDs 
are all independent of each other, and describe different aspects of the nucleon structure. However, it is natural to 
encounter relations among TMDs in simple models limiting to the valcncc-quark contribution and imposing SU(6) 
symmetry. The specific form of the relations can be traced back to the Melosh rotations which relate longitudinal and 
transverse nucleon polarization states in a Lorentz- invariant way. A similar situation occurs with the bag model [76] . 
In the diquark spectator model of Ref. [75] the relations (31) and (32) hold only for the separate scalar and axial 
contributions, while Eq. (33) is verified more generally for both u and d flavors. Since only two out of the four functions 
/ii giL, hi, h^rp are linearly independent, there are numerous relations among them. For example, subtracting (31) 
and (32) one gets a particularly interesting relation between pretzelosity, transversity and helicity distribution [76] 

5?(x,pt) - hUx,PT) = hf^'^''{x,pT). (34) 

This relation was recently discussed also in connection with the quark orbital angular momentum distribution [79] . In 
the version of the diquark spectator model of Ref. [78] the relation (34) is not supported in the axial-vector diquark 
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sector, but it remains valid for the scalar sector (sec also [79, 80]). Interestingly, in Ref. [81] the h^j,{x,pT) distribution 
was reconsidered also within a covariant parton model with the remarkable finding that the model satisfies the relation 
(34) without assuming SU(6) symmetry. 

The results in Eqs. (24)- (29) are applied in the following to a specific CQM taking the form of the momentum wave 
function from Ref. [102] 



^({x„p^J) = 2(2^)3 



1 



-,1/2 



Mq X 1X2X3 



N' 



(A/2+/32)7' 



(35) 



where ui^ is the free-quark energy and N' is a normalization factor such that J d[X]\ip{{xi}, {p±i})\^ ~ 1. In Eq. (35), 
the scale /?, the parameter 7 for the power-law behaviour, and the quark mass m are taken from Ref. [102], i.e. 
P = 0.607 GeV, 7 = 3.4 and m = 0.267 GeV. According to the analysis of Ref. [103] these values lead to a very good 
description of many baryonic properties. 

The results Eqs. (24)-(29) are general and can be applied to any CQM adopting the appropriate nucleon wave 
function. For example, we also considered the prediction in the hypercentral CQM model of Refs. [104, 105]. It has 
been observed that the description of nucleon properties using the model wave function either from [103] or from 
[104, 105] agree typically within (10-20) %, which might be considered as an indication of the accuracy of the CQM 
approach. In the following we shall assume that such is also the accuracy of the T-even TMDs from CQM [103]. The 
numerical results for T-even TMDs obtained in this way were discussed in detail in Ref. [77]. In order to compute 
T-odd TMDs it is necessary to go beyond the mere CQM scenario, and introduce gauge-boson degrees of freedom, 
which was beyond the scope of Ref. [77] and this work, where we concentrate on asymmetries due to T-even TMDs. 

In Fig. 2 we show the results for the integrals mpj, of the TMDs defined in Eq. (22), omitting the flavour dependence 
given by the SU(6) isospin factors iV" and P" in Eqs. (29). The solid curves correspond to the total results, obtained 
as the sum of the partial-wave contributions. The other curves show the contributions of the different orbital angular 
momentum components of the nucleon wave function. The unpolarized distribution /i, the hclicity distribution gi, 
and the transversity hi involve matrix elements which are all diagonal in the orbital angular momentum. In the 
plots of these functions, the dashed curves give the contribution from the S-wave component, and the dotted curves 
correspond to the P-wave contribution. The D-wave gives a negligible contribution to the /i and gi distributions 
(dashed-dotted curves), while it is absent in the case of hi. Although all these three functions are dominated by the S 
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FIG. 2: Parton distribution functions and transverse moments of TMDs as functions of x from the light-cone CQM [77]. In 
all panels the solid curves show the total results for the 'flavour-less' TMDs, i.e. the TMDs of definite flavour follow from 
multiplying by the spin-flavour factors A*'" or P°, see Eqs. (24)- (28). The other curves show how much the different angular 
momentum components of the nucleon wave function contribute to the total results: In the case of fi{x), gi{x), hi{x) the 
dashed, and dotted curves correspond to the contribution from the squares of the S- and P-wave components of the nucleon 
wave function, respectively. The D-wave contribution is absent in hi, while for /i and gi it is quite small and corresponds 
to the hardly-visible dashed-dotted curves. In the case of gf'^^\x), h\^l^{x) the dashed and dotted curves give the results 
from the S-P and P-D interference terms, respectively. In the case of h']l^^{x), the dashed curve is the result from the P-wave 
interference, and the dotted curve is due to the interference of S and D waves. 
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waves, they have a non-negligible contribution also from the P waves, with the largest (smallest) weight in the case of 

/i (pi). The functions gj'^^'' and hfj^^^ involve a transfer of orbital angular momentum by one unit between the initial 
and final nucleon state. In our model, they are simply related by Eq. (33). For these functions, the dashed curves in 
Fig. 2 show the contribution from the interference of S and P waves, and the dotted curves correspond to the results 
from the P- and D-wave interference term. The S-P interference term gives the largest contribution in the full x range, 
while the P- and D-wave interference term contributes at most by 20%. In the case of hf^^^ one has a mismatch of 
orbital angular momentum between the initial and final nucleon state equal to AL^ = 2. In the plot of this function 
in Fig. 2, the dashed curve gives the result for the interference of the = 1 and Lz = —1 components, and the dotted 
curve refers to the contribution from the interference of the the S- and D-wave components. Thanks to the interference 
with the S wave, we note that here the contribution from the D wave is amplified. Furthermore, at variance with 
the other distribution functions, the different partial-wave contributions do not have the same x-dependence, and for 
X > 0.6 the P waves are suppressed with respect to the S-D wave interference term. This peculiar behaviour makes 
the hf^^"^ function interesting, especially in the study of the interplay between the different partial wave components 
in the azimuthal spin asymmetries, as discussed in the Appendix. 

IV. COLLINEAR DOUBLE SPIN ASYMMETRIES Ai AND All 

Before discussing azimuthal asymmetries in SIDIS, we consider first the double spin asymmetry All and its inclusive 
analog Ai, Eqs. (13, 14). The study of these observables in the model framework is instructive, because in this case 
evolution equations (and fragmentation functions) are known and complications due to pT-dependence are avoided. 
This allows us to test the model under 'controlled conditions' in two respects. First, in which x-range and with what 
accuracy is the model applicable? Of course, the performance of the model could vary with observables. Nevertheless, 
this exercise will give us valuable insights in this respect. Second, how stable are the results under evolution? In this 
case we can compare exact results, with results obtained making assumptions on the evolution. The experience made 
here will be useful later, when dealing with azimuthal asymmetries whose evolution is practically not solved. 

A related key question emerging not only here but in any nonperturbative calculation concerns the scale at which 
the model results for the parton distributions hold. From the point of view of QCD where both quark and gluon 
degrees of freedom contribute, the role of the low-energy quark models is to provide initial conditions for the QCD 
evolution equations. Therefore, we assume the existence of a low scale Qq where glue and sea quark contributions are 
suppressed, and the dynamics inside the nucleon is described in terms of three valence (constituent) quarks confined 
by an effective long-range interaction. In fact, glue and sea quark degrees of freedom might be thought of at this low 
scale to be contained in the structure of the constituent quarks, which are massive objects. The actual value of Qq is 
fixed evolving back unpolarized data, until the valence distribution matches the condition that the second moment, 
i.e. the momentum fraction carried by the valence quarks, is equal to one [106]. 

Using LO evolution equations, we find Qq = 0.079 GeV^ [10'''] • Although there is no rigorous relation between 
the QCD quarks and the constituent quarks, and a more fundamental description of the transition from soft to hard 
regimes would be very helpful, this strategy reflects the present state of the art for quark model calculations [108-110], 
and has been validated with a fair comparison to experiments [106, 108]. 

Fig. 3a shows the inclusive double spin asymmetry Ai in DIS off proton, Eq. (14). The two theoretical curves are 
obtained using gf{x) and /f (x) from the light-cone CQM [77]. In one case both distribution functions are LO-evolved 
from the low scale of the model Qq to = 3.0 GeV^ (solid curve) using the evolution codes of Refs. [Ill, 112], and 
in the other case both distribution functions arc taken at the low scale of the model (dashed curve) . The results differ 
moderately at a; > 0.1 reflecting the weak scale dependence of Ai [113]. 

As can be seen in Fig. 3a, the description of data from the E143, EMC and SMC experiments [114-116] is reasonable. 
For X > 0.15 the model describes the Ai data within an accuracy of about 30%. The description improves in the 
valence- a; region of x > 0.2 though the accuracy of the data at large x is not sufficient to draw definite conclusions. 
Since the model contains no antiquark- and gluon-degrees of freedom, it is not surprising to observe that it does not 
work at small-x. As an intermediate summary, it can be said that the results are weakly scale dependent, and the 
model well catches the main features of the observable Ai in its range of applicability, namely in the valence-x region. 

Since in the following we will deal with SIDIS, we repeat the exercise with the double spin asymmetry All, Eq. (13). 
Figs. 3b, c show All in DIS production of charged hadrons from proton. The theoretical curves are obtained using 
gi{x) and /f (x) from the light-cone CQM [77], once LO-evolvcd to = 2.5 GcV^ (solid curve), and once left at the 
initial scale of the model (dashed curve). For the fragmentation function D^{z) we use in both cases, and throughout 
this work, the LO parametrization [117] at = 2.5 GeV^. Again we observe a weak scale dependence, and a good 
description of data in the valence x-region where the model describes the data within an accuracy of 0(20%). Thus, 
in the SIDIS case we make comparably positive experience as in the inclusive case. Notice that the result with gi{x) 




FIG. 3: The inclusive (a) and semi-inclusive (b, c) double spin asymmetries, Ai and All, defined in Eqs. (13, 14), in DIS off 
proton as functions of x. The theoretical curves are obtained with gi{x) and fi{x) from the light-cone CQM [77] as follows: 
both functions LO-evolved to the (Q^) of the experiments (solid curves), and both at the low scale of the model (dashed curves). 
In (b, c) we use always the parametrization [117] for at — 2.5 GeV^. The data in (a) are from Refs. [114-116], in (b,c) 
are from SMC (open circles) [118] and HERMES (black squares) [119]. 



and fi{x) taken at the low scale is, strictly speaking, not the consistent result for All at such a low scale because 
we use the parametrization for Z?" at = 2.5 GeV^. 

We remark that we could have tried to describe the double spin asymmetries Ai and All with g^x) from the 
model, and fi{x) from a parametrization, for example [120]. Such an approach would correspond to the strategy 
to use the model only as input for the part which is responsible for the spin effect, and to use for the well known 
denominator of the spin asymmetry standard parametrizations, which has the advantage that the model uncertainty 
is only in the numerator. In the case of Ai and All, however, such an approach yields a bad description of the data. 
This can be traced back to the fact that the /f (x) from the model [77] and from parametrizations [120] have different 
large- a; behaviour. Interestingly, it happens to be the case also in the case of gf{x) from the model [77] and from 
parametrizations [121], such that the uncertainties partly cancel in the ratio, leading to a better description of the 
data. It is important to stress that here we deal with chiral-even functions, where antiquark and giuon degrees of 
freedom are of importance. In the case of chiral-odd TMDs the situation is different, and a different approach could 
be more successful. We will come back to this point later on. 

The above discussion allows to assign a 'typical accuracy' to the approach. In this context it is of interest to make 
the following observation. The present version of the model uses SU(6) symmetry, such that gi{x) = —Agf{x) at the 
low scale, and similarly for other polarized distribution functions or TMDs. Due to isospin symmetry the structure 
functions of the neutron are related to those of the proton by interchanging u and d flavour. Therefore, at the low 
scale Ai cc | -I- | gf^^ = | 5i + | ffi = 0- At higher scales A" ^ due to evolution but the effect remains small. 
Rather than claiming A" « 0, it is more meaningful to state that SU(6) predicts A" to be small compared to, say, 
the Ai of proton. Thus, for the SU(6) symmetry concept to be a useful tool, we expect for the ratio 

A" 

«1. (36) 

Experimentally A" (extracted by subtracting deuteron and proton data, or from ^He data, modulo nuclear corrections) 
is found clearly non-zero. However, in the valence-x region the SU(6) expectation is supported by data [122-127]: the 
ratio (36) is of the order of magnitude 20 % — which is indeed a 'zero' within the model accuracy. 

Notice that in SIDIS the SU(6) symmetric TMDs are weighted with fragmentation functions, such that in general 
azimuthal asymmetries from the neutron are non-zero already at the low scale (with the exception of tt" into which 
u and d quarks fragment with equal strength). Nevertheless, also in SIDIS the results for a neutron target are 
highly sensitive to SU(6) breaking effects, and we refrain from showing them here. An adequate description of spin 
asymmetries from a neutron target requires to account systematically for possible SU(6)-breaking effects. This is 
similar to what one observes in the case of the electric form factor of the neutron, where S'-wave components in the 
nucleon wave function were found essential to reproduce the experimental data [85, 128]. Results for the TMDs with 
such SU(6) breaking terms will be discussed elsewhere. 

Let us draw conclusions from the study presented in this Section. The double spin asymmetries in inclusive DIS, Ai, 
and SIDIS, All, sue weakly scale-dependent. The model describes the data on these observables within an accuracy 
of 20 — 30% in the valence-x region. This suggests that the approach could also be useful for studies of azimuthal 
spin asymmetries. In fact, at this point it is worth to stress that (i) azimuthal phenomena are expected to yield sizable 
effects especially in the valence- a; region, (ii) Erst data on azimuthal asymmetries often have uncertainties comparable 
to the observed model accuracy, (iii) in proposals for future experiments predictions of new effects of an accuracy of 
0(30%) are useful enough. 
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V. Pt -DEPENDENCE, GAUSS ANSATZ, APPLICABILITY OF THE MODEL 

When dealing with azimuthal asymmetries in SIDIS it is very convenient to use the Gaussian model for the distribu- 
tion of transverse parton momenta, see Sec. II. If one assumes the Gaussian Ansatz (15) for fiix^px) and D^{z, Kt), 
then a good description of the SIDIS data (more precisely, mean values for {Pi^±{z)) not corrected for acceptance 
effects) from HERMES [33] is obtained with the following parameters [58] 

(Pt(/i)) = 0.33 GeV^ , {Kt{Di)) ^ 0.16 GeV^ . (37) 

Numerically very similar results were obtained in [56] from a study of EMC data [30] on the Cahn effect [5] . 

The pT-dependence of TMDs in the model [77] is definitely not of Gaussian form. However, the essential question 
is: Can it be reasonably approximated by a Gaussian form? 

In order to discuss that let us make the following two exercises. First, wc define the mean transverse momenta 
(n ~ 1) and the mean square transverse momenta {n ~ 2) in the TMD j{x,pT) as follows 

'■^ J da; J d^pT j{x,Pt) 

In Table I we show results for these quantities for T-even twist-2 TMDs from [77]. 

In order to see to which extent the results for the pT-dependence of TMDs from [77] can be approximated by a 
Gaussian behaviour, we remind that in the Gaussian model the following relation holds 

(39) 

TT 

In Table I wc show also the results for the ratio that would be unity for a Gaussian pr-distribution. Remarkably, 

the model results for this ratio from [77] deviate from unity by not more than 10 %. Of course, although the Gaussian 
model relation (39) works within 10 %, it does not necessarily imply that the pT-dependence in the model is Gaussian 
within such an accuracy. We make therefore the following second exercise. 

We ask the question: what is the difference between computing in the model [77] an observable using exact model 
PT-dependence of TMDs and computing it by approximating the true pT-dependence by a Gaussian? We can rephrase 
this question also as follows: when integrating out the transverse momenta of produced hadrons and focusing, for 
example, on the x-dependence of azimuthal asymmetries, the pT-model-dependence is weakened, but to what extent? 

In order to answer that question we choose the double spin asymmetry ^™^('^''~'^^) aj^j ^ge f]^^ model [83] for 
Di{z,Kt), which also refers to a low scale. (We stress that the results presented here are to be considered as an 
exploratory study of py-model effects. Our predictions for this asymmetry will be given in the subsequent Section.) 

In Fig. 4 we see the results for obtained as follows. The solid curve shows the result from solving 

numerically the convolution integral in (6) with g^rp{x,pT) from [77] and D1{z, Kt) from [83]. The dotted curve 

shows the result for the asymmetry obtained from the Gaussian model, Eq. (17), using g^rj^^'^"' (x) from [77] and D1{z) 
from [83], assuming the Gaussian Ansatz (15) for these functions, and assigning the Gaussian widths according to 
Eq. (38). 



inGeV inGeV^ 7r(p2> (Pt(/i)> 



fi 0.239 0.080 0.909 1 

ffi 0.206 0.059 0.916 0.74 

hi 0.210 0.063 0.891 0.79 

gir 0.206 0.059 0.916 0.74 

0.206 0.059 0.916 0.74 

hir 0.190 0.050 0.919 0.63 



TABLE L The mean transverse momenta and the mean square transverse momenta of T-even TMDs, as defined in Eq. (38), 
from the hght-cone CQM [77]. If the transverse momenta in the TMDs were Gaussian, then the result for the ratio in the 
fourtli column would be unity, see text. The last column shows the {ptU)) in units of (Pt(/i)). 
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The different results agree within an accuracy of 20%, see 
Fig. 4. Such an uncertainty is 'within the model accuracies' 
of Refs. [77, 83]. Thus we conclude that the true transverse- 
momentum dependence in the models [77, 83] can be approxi- 
mated by the Gaussian Ansatz with a satisfactory precision for 
practical purposes. 

Next we address the question how to use consistently the 
model predictions [77] for phenomenology — in view of the fact 
that they refer to a very low hadronic scale. The fact that 
of in that model is smaller compared to what is required by 
phenomenology, Eq. (37), is perfectly reasonable. Sudakov ef- 
fects make the p^-distributions broader, i.e. (p^) larger, when 
evolving to larger (experimentally relevant) scales. 

This j>T-broadening is expected to be independent of the quark 
polarization, in first approximation. Thus, what we can use for 
phenomenology are the model results for (p^) in units of the 
mean square transverse momenta of /i, see last column in Ta- 
ble I, and take the 'unit' (Prifi)) from phenomenology, Eq. (37). 

On the basis of the considerations in this and in the previous 
Section we are in the position to establish our strategy to treat 
azimuthal asymmetries in the following. Let us summarize. 

• We will mainly focus on the a--dcpendence of the asymmetries, especially in the valence-x region (sec Sec. IV). 

• We will assume the Gaussian model, which is a reasonable approximation (this Section, see above). 

• When information on a specific Gaussian width of a polarized TMD is needed, we will use the model prediction 
for the corresponding ratio (see last column in Table I), and the value from Eq. (37) for the width of . 

• We will not discuss the z-dependence of the azimuthal asymmetries, because here integrals over the x-dependence 
enter which extend, depending on the experiment, to low- a; regions where the model is not applicable.^ 

• Similar warnings apply to the P^^-dependence of the asymmetries. We shall therefore address this point with 
particular care, see Sec. X below. 



0.1 



0.08 




0.02 



0.2 0.4 0.6 X 

FIG. 4: ^^°^('*'>^'^s) ^+ pi-oduction off proton, as 
function of X. Solid curve; exact result obtained using 
girixypr) from [77] and D\{z,Kt) from [83]. Dashed 
curve: an approximation obtained using the inte- 
grated functions (s:), D\{z) from [77, 83] and 'sim- 
ulating' their pT-dependence by means of the Gaus- 
sian Ansatz, as described in the text. 



VI. THE DOUBLE-SPIN ASYMMETRY ^=°^(*h-*s) 



We start the discussion of azimuthal asymmetries with the double-spin asymmetry A^"^^^'^'' — F^^^'^'" '^^^ /Fun 
which is proportional to g^rji 75°, see Eq. (6). Assuming the Gaussian Ansatz, which gives a good approximation, 

see Sec. V, we have to model the prefactor B'q in Eq. (17). For that let us rewrite that factor as 



-1/2 



For the first ratio in the curly brackets we use the model prediction from the last column in Table I. For the second 
ratio in the curly brackets we use the numbers from Eq. (37). 

Being interested in the x-dependence of the asymmetry, we further integrate over z 



,cos(0,^0s),^^ _ Y.a(iixg^^^^''{x) {B'„D1) 



where (. . .) denotes the average over z within the respective experimental cuts. Here and in the following, we will 
consider the range 0.2 < z < 0.7 corresponding to the typical kinematics of HERMES. There is little difference if one 



We recall that the numerators and denominators of the asymmetries (5)— (10) are actually weighted by 1/Q* oc which strongly 

emphasizes the role of the small- a; region, whose description is beyond the range of applicability of the model. 
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FIG. 5: The double-spin asymmetry A'^^p in DIS production of pious, as function of x, obtained using (jjy (s) and fi{x) 

from the light-cone CQM [77] in the following way: both functions are taken at the low scale of the model (dashed curves), and 
both are LO-evolved to = 2.5 GeV^ (solid curves). Hereby the scale dependence of g[}jl'^{x) is 'simulated' using the gi{x) 
evolution pattern, see text. The data points are preliminary COMPASS data for charged hadron production off deuteron [45]. 



uses the COMPASS cuts 0.2 < z < 1, since the resulting (z) is similar. At JLab typically higher (z) are reached. For 
the present observable, however, this has little impact. 

The results for -^^ DIS-production of pions off different targets are shown in Fig. 5. For gj^^^°(x) and 

fi{x) we take the results from the model [77], and consider two options. First, we take both functions at the low 
scale of the model (dashed curves). Second, we consider both curves LO evolved to ~ 2.5 GeV^ (solid curves). 

Hereby, we use for g^^^'^°'{x) the evolution equations for .gf (x). This is admittedly not the correct evolution pattern. 
However, this is the evolution pattern of a chiral-cven polarized function, and the purpose of presenting it here is to 
shed some light on the possible size of evolution effects. 

Our crude estimate of evolution effects indicates, that the predictions for the asymmetries are presumably robust 
concerning scale dependence. The proton asymmetries reach 4 % in the valence-x region, which could be measured 
— especially at JLab. The deuteron asymmetries are somewhat smaller. For a deuteron target there also exist 
preliminary data from the 2002-2004 run of the COMPASS experiment [45]. As can be seen in Fig. 5, our results are 
compatible with these preliminary data. 

Estimates for A'^^^^'^'' were made also in [66] on the basis of the approximation 

,-«''(x,g^)"«^xr^5?(y,Q^) (42) 
Jx y 

using the parametrization [121] for gi{x). The approximation is 'justified' in QCD upon the neglect of pure twist- 
3 (quark-gluon) correlators and current quark mass terms [67-69]. This is analog to the Wandzura-Wilczek (WW) 
approximation for the twist-3 parton distribution function g^{x) [129-131] — hence the label 'WW' in (42). The WW- 
approximation for g!^{x) is supported experimentally within the error bars of the present data [132-134]. Whether 
the WW-type approximation (42) is supported by data equally well remains to be seen. 

In the light-cone CQM [77] the 'WW-type approximation', Eq. (42), is supported in the valence x-region with good 
accuracy. Furthermore our results support the findings of [66] also numerically. Taking into account the different 
kinematical cuts applied in the calculation of [66] , we obtain asymmetries of similar size, with a more flat x-dependence. 
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VII. THE SINGLE-SPIN ASYMMETRY 

Next we focus on the azimutlial SSA /^^'^'^^'^^^ = F^^'^'^^'^^'^ / Fuu due to transversity and the CoUins function. 
In the Gauss Ansatz (15) the structure function in the numerator of this SSA is given by the expression in Eq. (18). 

For the CoUins function, more precisely for {BiH^''^^^^'^) equal to (2i?Gauss^^i^'^^^'''') in the notation of [61], we use 
the results extracted in [61] from the (preliminary) HERMES data [37]. Although meanwhile new data are available 
[43, 44, 46, 47] the results on H^^ from [61] are still in excellent agreement with updated extractions [63]. 

When dealing with asymmetries due to chiral-odd TMDs, in our opinion a different approach is more appropriate 
as compared to the case of asymmetries due to chiral-even TMDs. Let us explain this point in more detail. 

When describing asymmetries due to chiral-even functions in the previous Sections, we used model input for both, 
the numerator and the denominator of the asymmetries. In the model gluon (and sea quark) degrees of freedom are 
absent at the low scale, and generated by evolution at higher scales. Admittedly, in this way one cannot accurately 
describe absolute DIS cross section data. For that non-zero unpolarized gluon and sea quark distributions are needed 
already at low input scales [120] (though the model scale is lower than the initial scale of the parametrizations [120]). 
This 'shortcoming' of the model, however, affects similarly the numerator and the denominator of asymmetries due 
to chiral-even TMDs. Indeed, we observed that these model uncertainties partly cancel in the ratio — resulting in a 
useful description of (SI)DIS data on asymmetries in the valence-x region, see Sec. IV. 

Can we expect a similarly good description of asymmetries, which are due to chiral-odd TMDs, using this strategy? 
The answer is no, in our opinion. Transversity has no gluon counterpart, in contrast with fi{x). The absence of 
gluon degrees of freedom in an approach constitutes therefore a 'lesser shortcoming' for hi{x) than for /f (x). So one 
expects intuitively that in quark models transversity and other chiral-odd TMDs could be modeled more reliably than 
chiral-even ones, though it is not clear how to put this expectation on a firm field theoretical basis. 

Nevertheless, these considerations suggest to adopt the following strategy for the description of the Collins SSA, 
namely to use hi{x) from the model LO-evolved [55, 135] to the experimental scale in the numerator of the SSA, 
and fi{x) from a parametrization, e.g. [120], taken at the corresponding scale. In this way, the model uncertainty 
is limited to the numerator of the SSA only, while the denominator is described exactly. We indeed observe that 
the above-described strategy yields by far the best results in the case of the Collins SSA, see Fig. 6. Any other 
options, such as hf{x) and /f (x) from the model at the low scale or hf{x) and fi{x) from the model LO-evolved, 
gave unsatisfactory results. 

Let us discuss in some more detail the results for A'j^'^'^'^'^'^^'' in Fig. 6, where for sake of clarity we refrain from 
showing the error bands due to the statistical and systematic uncertainties of the extracted Collins function [61]. 
Figs. 6a, b show the results for charged pion production from a proton target in comparison to the preliminary 
HERMES data [37]. (It is consistent to compare to these data, because the information on H^^ [61] was extracted 
from those data.) The model results ideally describe these data — including the small- x- region, sec Figs. 6a, b. This 
is in line with the favourable comparison between our model predictions [55, 101] and the phenomenological extraction 
of the transversity and tensor charges in Ref. [62, 63]. 

In Figs. 6c, d we compare our results for A^^"^^ = —A^'^'^'^'^'^^'^ (since if)c — 4>h + 4>s + ^) for charged pion 
production from a deuterium target to the COMPASS data [46], which extend down to much lower values of x. Our 
results are compatible with the data also in this case, including again the small- .t region. 

On the basis of the presently available information on the Collins function extracted from SIDIS and e+e^ data, 
we would predict tt" SSAs compatible with zero within the uncertainties of the extractions [61-63]. However, this is 
a prediction due to our present understanding of the Collins effect, rather than due to the model for TMDs. For this 



A^'*h+*s' (a) A™'*h-^*s» (b) A'™^*c' (c) A™i*c) (d) 




FIG. 6: The single-spin asymmetry Aj)'^''*''^'^^^ = —A^^^'^ in DIS production of charged pions off proton and deuteron 
targets, as function of x. The theoretical curves are obtained on the basis of the light-cone CQM predictions for hi{x,Q^) 
from Ref. [55, 77], see text. The (preliminary) proton target data are from HERMES [37], the deuteron target data are from 
COMPASS [46]. 
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reason, here and in the following two Sections where we discuss further SSAs due to Collins effect, we refrain from 
showing results for neutral pion production. 



VIII. THE SINGLE-SPIN ASYMMETRY A^^' 



In this Section we discuss the azimuthal SSA A^'^^^'^''^ = F^^^'^'^''^ / Fuu due to h^f; and the Colhns function. In the 
Gauss Ansatz (15) the structure function in the numerator of this SSA is given by the expression in Eq. (20). Thus, 
in order to describe this SSA we need {Bl^H^'^^^'^^"') which we estimate on the basis of the Collins function extractions 
[56; 58] precisely as described in Ref. [67]. 

The problem we face in the context of concerns the question how to evolve correctly h^^''°'{x) from the low 

initial scale of the model to the relevant experimental scale. In contrast with transversity, exact evolution equations 
are not available in this case. 

In our study of A!^^'^^'^'^^'^ we learned that other strategies, such as leaving transversity at the low scale of the 
model (and taking fi{x) in the denominator from the model or from parametrization, at the low scale or evolved) 
resulted in unfavourable descriptions of data, and we were able to understand qualitatively why. Of course, this is 
here a different observable. But the experience with the Collins SSA does not encourage any other strategy than that 
adopted in that case, in Sec. VII, namely to evolve the chiral-odd TMD from the model, and use parametrizations 
for the denominator of the SSA. 

Not being able to evolve hfj^^^°'{x) correctly we use instead the /i°(a;)-evolution pattern to evolve it 'approximately'. 
Since both functions are chiral-odd, the simulation of evolution effects in this way can be expected to be more 
promising than using any other evolution pattern. 

In Fig. 7 we compare the results obtained in this way to the HERMES data from proton and deuteron targets 
[31, 33]. We observe that our estimates are well compatible with the data — including again the small- a; region. For 
comparison, in Appendix and Ref. [136] we also show the results obtained with the same ingredients as in Fig. 7 but 
without approximate evolution of h1^{x). This approach yields a somewhat larger SSA, especially at large x, but it 
is similarly compatible with the data. 

In Ref. [67] predictions for the ^^^(^^fi) were made on the basis of the WW-type approximation 

h-i^^\x)^^-x^ f%K[y), (43) 

J X y 

and model predictions for transversity from [137]. Eq. (43) is analog to the approximation (42), i.e. it also arises 
when certain quark-gluon correlator and current quark-mass terms arc neglected. Interestingly, the light-cone CQM 
supports the approximation (43) within a reasonable accuracy [77]. Also the numerical results for the SSA obtained 
here and in [67] agree well qualitatively. 

It is, of course, an important question how to quantify the theoretical uncertainty we introduced in our study 
by employing the incorrect evolution pattern for h^^^°'{x). Until exact evolution equations for this TMD will be 
available, this question cannot be answered exactly. However, one may suspect that the uncertainties due to evolution 
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FIG. 7: The single-spin asymmetry A^^^'^'^'^^ in DIS production of charged pions off proton and deuteron targets, as function 
of X. The theoretical curves are obtained by evolving the light-cone CQM predictions for h^^^'^" of Ref . [77] to = 2.5 GeV^ 
using the hi evolution pattern, see text. The data points are from HERMES [31, 33]. The inner error bars are the statistical 
errors, the outer error bars are the systematic errors. 



14 
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FIG. 8: The single-spin asymmetry A^'j; in DIS production of charged pions off proton and deuteron targets, as function 

of X. The theoretical curves are obtained by evolving the light-cone CQM predictions for /if^l^'" of Ref . [77] to = 2.5 GeV^ 
using the evolution pattern, see text. The preliminary COMPASS data are from Ref. [45]. 

are less dominating than other uncertainties within the model. The current HERMES data do not contradict this 
expectation, see Fig. 7. We remark that there are also preliminary CLAS data [40]. Our approach is compatible with 
the results for tt^ and 7r° but cannot explain the trend of the tt~ SSA, similarly to Ref. [67]. The situation will be 
fm-ther clarified in future experiments at JLab [138, 139], and COMPASS. 

IX. THE SINGLE-SPIN ASYMMETRY 74^^'^*""*''^ 

Finally we study the azimuthal SSA A™^^^'^'^"'^^' = i^^^^'^'^+'^^V-F'c/c/ due to pretzelosity and the Collins function. 
In the Gauss Ansatz (15) the structure function in the numerator of this SSA is given by the expression in Eq. (21). 

The factor {B^H^^^^^^"') we evaluate exactly as done in Ref. [76]. 

Also in the context of the asymmetry A'^j^'^^'^'^ ^^"^ we face the question how to evolve hf^^^°'{x) from the low initial 
scale of the model to the relevant experimental scale. Exact evolution equations are not available in this case, either. 
We follow here the approach developed in the previous Section, and 'simulate' the evolution of hf^^'''^{x) by evolving 
it according to the transversity-evolution pattern. Again, since pretzelosity and transversity are chiral odd, this way 
might be a useful estimate of evolution effects. 

The results obtained in this way are shown in Fig. 8. We find the pretzelosity SSA rather small, about one percent 
in the case of charged pions from a proton target, see Figs. 8a, b. This makes it the probably most challenging 
asymmetry to be measured. The deuteron SSAs are somewhat smaller, see Figs. 8c, d where we show for comparison 
the preliminary COMPASS data presented in Ref. [45]. Our results are compatible with the data, and explain why 
the effect was found consistent with zero within error bars at COMPASS. The error bars of the preliminary data [45] 
simply do not allow to resolve an asymmetry smaller than one percent. 

In Ref. [76] estimates for the asymmetry A^^'^'^''"'^^'' were presented on the basis of the positivity bound 

l/i;f^^^°(2;)l < i(/f(a;) - g'^{x)) [48], using the parametrizations [120] for f^{x), gf{x). The results of the light- 
cone CQM for pretzelosity (as well as other TMDs), of course, respect positivity bounds [77], and the transverse 
moment of pretzelosity at the low scale of the model is not that small, see Fig. 2. But after evolution (with the 
transversity evolution pattern) to a scale of = 2.5 GeV'^, it is much smaller than its bound constructed from pa- 
rameterizations for fi{x), gi{x) at = 2.5 GeV^. Therefore, our estimates of the pretzelosity SSA are significantly 
smaller than the maximal effect allowed by positivity requirements [76] . 

Of course, we do not know to which extent our approach to estimate the hf^^\x) evolution effects is really realistic. 
For comparison, in Appendix and Ref. [136] we also make predictions neglecting the evolution of pretzelosity. In this 
way the results for the SSA are more sizable, and the effects are larger especially in the region of intermediate and 
large x. The planned experiment at JLab will allow us to discriminate among the different predictions [140]. 

X. Ph_L -DEPENDENCE OF SPIN ASYMMETRIES 

As discussed in Sec. V, care is required in order to use the model results for the pT^-dependence of TMDs for 
phenomenological applications. In this Section we shall exemplify how this can be done with a study of the Ph±- 
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dependence of the double spin asymmetry All. In principle, we could discuss also other asymmetries, but All has 
the advantage that its x and z dependence is rather well known — so we do not need the model input for that, 
and can focus on P/j^-dependence which is the only new concept in this case. Would we discuss other (azimuthal) 
spin asymmetries, we would need to use the model input also for the x-dependence of the novel TMD, and face the 
problems of how to evolve TMDs, make a meaningful estimate of Sudakov effects, and deal with the small-x region 
(see footnote 1). When dealing with All in the way described below, we avoid these problems. 

Before discussing the P/iX-dcpendcncc of All in our approach, let us remark that ideally a study of px-effects 
should start with absolute cross section data on the production of hadrons in unpolarizcd DIS. Such data arc difficult 
to produce, and experimentally it is preferable to study the Pft,_L-dependence of asymmetries, since detector acceptance 
effects in the numerator and denominator of the asymmetries (largely) cancel. Therefore, so far information from 
SIDIS on pr-dependence of the unpolarized parton distribution and fragmentation functions, /f and Df, has been 
obtained only indirectly, see Sec. V. It would be desirable to improve this situation. Apart from the absolute cross 
section proportional to Fuu, the next 'simplest' observable to learn about pT-effects is probably All = Fll/Fuu- 

In Sec. V we learned that the Gauss Ansatz is supported within the model with reasonable accuracy. This justifies 
to make explicit use of it, also in this case. If we assume this Ansatz, then 
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If we assume that the widths are flavour and 
asymmetry is given by 



or z-independent, then the P/,^-dependence of the double spin 
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where (All) denotes the spin asymmetry averaged over x and z, which is known in the experiment with precision. 
(Now in Eq. (46) it is implied that "z^" in (44), (45) is replaced by (z^). This is an approximation, and the treatment 
could be improved, but we refrain from this in our illustrative study for sake of clarity.) 
We remark that positivity, i.e. All < 1 V a; 



and Ph±, dictates 



{plih)) 



< 



(47) 



1.5 - 



This implies that (in the Gauss Ansatz) the 
widths of helicity and the unpolarized distribu- 
tion could be equal, if and only if the equality 
ff{x) = \gi{x)\ were true. 

Now, let us discuss how to use the model re- 
sults in order to predict the P;ij^-dcpendcncc of 
All- From Table I we know that 



(pUh)) 



= 0.74 



(48) 
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If we take this ratio for granted, and assume Q 0,5 

for (p|(/i)) the result from [58], Eq. (37), then 
we obtain for a{Phi_) = ALL{Ph±)/ (All) the 
results shown in Fig. 9. We include in Fig. 9 also 
predictions based on using the results for (K^), 
(Pt(/i)) fr'om [56]. We observe a rather stable 
prediction which depends little on the choice of 
parameters [58] vs. [56]. The prediction in Fig. 9 depends more strongly on the model prediction (48). 

This result is (in our approximations) the same for any target and produced hadron. In fact, in the SU(6) symmetric 
light-cone CQM of Ref. [77] the widths as defined in Eq. (38) are always fiavour independent. But wc recall, that the 
entire Gauss Ansatz is in the light of the results of Ref. [77] merely an approximation. 



Pfjj_ in GeV 

FIG. 9: a{Ph.±) = ALL{Ph±)/ {All) vs. Ph± in SIDIS, for experiments 
with (z^) = 0.16. The results are obtained using the prediction (48) 
from the model [77] , and the Gauss model parameters in Eq. (37) from 
Ref. [58] (solid curve) or the corresponding parameters from Ref. [56] 
(dashed curve). 
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It will be instructive to learn to which extent our predictions will be confirmed by experiment. As mentioned, 
we could similarly discuss predictions from the model for azimuthal asymmetries, too. But those predictions would 
presumably have larger theoretical uncertainties, such that we shall content ourselves here with the study of the 
Pft_L-dependence oi All- 

XI. CONCLUSIONS 

In this work we have studied all leading-twist azimuthal spin asymmetries in SIDIS due to T-even TMDs on the 
basis of predictions within one and the same model, i.e. the light-cone CQM of Ref. [77]. 

By studying first the well known double spin asymmetries Ai in DIS and All in SIDIS, we demonstrated that 
the approach is capable of describing the data on these asymmetries in the valence- a; region with an accuracy of 
0(20 — 30) %. The comparison with results from other constituent models has shown this to be a typical accuracy to 
which the constituent quark model scenario can be expected to work. 

We paid particular attention to the question, how to apply the model results for TMDs obtained at a very low 
hadronic scale to the description of data referring to high scales of typically several GeV^. We made a test for the 
double spin asymmetries Ai in DIS and All SIDIS where the evolution equations involving the parton density 
fi{x) and the helicity distribution gi{x) are exactly known. In these cases we have been able to demonstrate the 
stability under evolution of our results in the valence- a; region. 

For TMDs entering the description of azimuthal asymmetries, however, not only the evolution with renormalization 
scale has to be taken into account, but also Sudakov effects which broaden the pT-distribution of the TMDs. We tack- 
led this issue in two steps. First, we observed that the light-cone CQM [77] do not show a Gaussian pT-dependence. 
Nevertheless, their pT-dependence entering the azimuthal asymmetries is integrated over in certain convolution inte- 
grals, so that we have found that within the accuracy of our approach the effect of the true pr-dependence in the 
TMDs can be approximated by a Gaussian dependence. We have therefore explicitly employed the Gaussian Ansatz, 
which allows to express azimuthal asymmetries in terms of parton distribution functions or transverse moments of 
TMDs. In the second step, we used evolution equations to evolve the respective parton distribution functions or 
transverse moments of TMDs to the experimental scales. 

We have been able to do this exactly, strictly speaking, only in the case of transversity h1{x). In the other cases, 
the evolution of the transverse moments of TMDs was estimated by employing those evolution equations, which seem 
most promising to be able to simulate the correct evolution, which is presently not available. For example, we evolved 
9iT^'^{^) by means of the evolution pattern of the (also chiral-even) gi{x), while for hf^^^°'{x), hf^^^°'{x) we used the 
evolution pattern of the chiral-odd h1{x). The theoretical uncertainties due to these approximate treatment of the 
scale dependence is presumably not larger than the accuracy of the model. 

Among the leading-twist azimuthal spin asymmetries due to T-even TMDs, the Collins SSA is the only 

non-zero one within the present day error bars. We observe a very good agreement of our results for the x-dependence 
of this SSA with the HERMES proton [35, 37], as weU as with the COMPASS deuteron target data [36, 39]. 

The presently available final data on A'^"^^^"^"^ [31, 33] or preliminary data on ^sMi^^~<t>s) j^gj 

show results compatible with zero within error bars. Our results are compatible with these first or preliminary data. 
In future, our predictions of these azimuthal spin asymmetries could be tested by more precise data — especially from 
COMPASS and JLab. 

In an exploratory study of the double spin asymmetry All we have shown how model results for TMDs obtained at 
very low scale could be applied for studies of the P/i^-dependence of spin asymmetries. We have chosen this observable, 
because here the P;i^-dependence is the only new aspect, the x- and z-dependence being known experimentally with 
good precision. For that we explored again the fact that the light-cone CQM [77] supports the Gaussian Ansatz 
for TMDs within a reasonable accuracy, and used as model input only the prediction for the ratio of the mean 
transverse- momentum squares of gl and /f . This ratio is expected to be little affected by Sudakov effects in a first 
approximation. We made predictions for ALL{Phi.) which could be tested soon, for example, at JLab [140]. 

The advantage of our study is that the same model input has been used to describe all leading-twist spin asymmetries 
due to T-even TMDs. Wherever the data allow to draw definite conclusions, we observed a good agreement with 
the experiment in the range of applicability of the approach. It remains to be seen whether also our predictions for 
the other azimuthal spin asymmetries will be similarly confirmed by future data. If so, our approach will provide 
interesting insights in the spin and orbital angular momentum structure of the nucleon, which — though being model 
dependent — are of interest by themselves, as it is exposed in the Appendix. 
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APPENDIX A: ANGULAR MOMENTUM DECOMPOSITION OF SPIN ASYMMETRIES 

111 this Appendix we discuss the contribution from the different angular momentum components of the nucleoli 
wave function to the spin asymmetries. To this aim, we calculate the numerator of the asymmetries using the results 
of the TMDs at the hadronic scale of the model, and separating them into partial wave contributions according to 
the decomposition shown in Fig. 2 of Sect. III. In order to discuss how this decomposition behaves under evolution to 
higher scale, one would need to know the evolution equations for the different angular momentum components of the 
nucleon wave function separately, but, to our knowledge, this problem has never been addressed so far and is beyond 
the scope of our work. Although this decomposition is model dependent and it is not possible to extract experimentally 
the absolute strength of the different partial waves, it is instructive to visualize how the angular momentum content 
of the TMDs affects the spin asymmetries. In particular, the combined analysis of different spin asymmetries can 
give insights about the relative strength of the different partial waves, and therefore can be useful in modeling the 
light-cone wave function of the nucleon. 

The results presented in this Appendix, corresponds to an alternative approach concerning the question how to 
use the model results referring to a low hadronic scale for phenomenology at experimentally relevant scales. Namely, 
here we use the model at the low scale only as input for the part which is responsible for the spin effects. For the 
well known denominator of the spin asymmetries we use standard parametrizations at the experimental scale. In 
this way, the model uncertainty is only in the numerator. The comparison of these results and those presented in 
Sees. VI, VIII, IX, where the attempts were made to approximate evolution effects of the TMDs, with forthcoming 
experimental data [138-140] may give us interesting information about the scale dependence of these observables. 

In Fig. 10, we show the results for the All asymmetry, obtained by using both the unpolarized distribution 
function /i and the helicity distribution function gi from the light-cone CQM at low scale. The total results are 
further split into the contributions to gi from the S- (dashed curve) and P-wave (dotted curve) components, while the 
D-wave contribution is not shown because it is negligible. These separate terms reflect the dominance of the S-wave 
component with respect to the P wave in gi, as already observed in Fig. 2 of Sect. III. Furthermore, the S-wave 
term is practically constant in the full x range, while the P-wavc contribution is slowly increasing at larger value of 
X, reaching a maximum of about 30% of the total result. 

The A™^^"^'' asymmetry shown in Fig. 11 is calculated with both the unpolarized distribution function /i and 

the helicity distribution function g[^r^^ from the light-cone CQM at low scale. Here we separate the contribution to 

pjy"'" from the interference of S and P waves (dashed curves) and from the interference of P and D waves (dotted 
curves). We see that the S- and P-wave interference term governs both the size and the shape in x of the total results, 
while the contribution of the P- and D-wave interference is rather small and constant in the full x range. 

We now pass to consider single spin asymmetries involving chiral-odd TMDs. In the following, we will use for /i 
the parametrization from Ref. [120] at = 2.5 GeV^, while for the hi, /i^T^' ^^'^ TMDs we will use the 

results from the light-cone CQM at low scale. The A^'^'^'^^'^^'^ asymmetry is shown in Fig. 12, with the separate 
contribution to hi from the S- (dashed curves) and P-wave (dotted curves) components. The contribution from the 
P waves is within 30% of the total results, and, at variance with the double spin asymmetries discussed above, the 
two partial-wave contributions have very similar x dependence, with a maximum at a; ~ 0.7. 

Ill Fig. 13 we show the results for the j\^^'^^^'^'^'> asymmetry, due to the to the h^j^ TMD. Since in our model 
h^j^ = —gix, the relative strength of the contributions from the S- and P-wave interference (dashed curves) and the 
P and D-wave interference (dotted curves) is the same as for the corresponding contributions in A™^'-'^''"'^^-'. 

Finally, in Fig. 14, we show the results for the asymmetry, separating the contributions to /ij^^ from 

the interference of P waves (dashed curves) and S-D waves (dotted curves). This is the only case where we can 
exploit the interference with the large S-wave contribution to amplify the effects due to the small D wave. The two 
interference terms have a quite different shape as function of x: in the case of the P-wavc interference, we have an 
oscillating behaviour, with a sign change at a; ~ 0.7, while the S-D wave interference term is similar to a bell-shaped 
curve with the maximum at a; ~ 0.7. The sum of these two contributions gives a total result which is peaked at 
X ~ 0.4. At larger x, the S-D wave interference term gives the main contribution, while at smaller x the P wave and 
the S-D wave interference terms contribute with the same strength. 
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FIG. 10; The double-spin asymmetry All in DIS production of charged pions off proton and deuteron targets, as function of x. 
The results are obtained using gf^j^^^'^lx) and fi{x) from the light-cone CQM [77] at the hadronic scale, and decomposing (/j'^l^'" 
into different partial wave contributions: the dashed curves correspond to the contribution from S waves, the dotted curves are 
the results for the P-wave contribution, and the solid curves are the total results, sum of the S- and P-wave contributions. 
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FIG. 11: The double-spin asymmetry yl™^^''^'' "^^^ in DIS production of pions off proton and deuteron targets, as function of x. 



The results are obtained using (7j^^-'"(a;) and fi{x) from the light-cone CQM [77] at the hadronic scale, and decomposing g^j, 
into different partial wave contributions: the dashed curves correspond to the contribution from the S- and P-wave interference, 
the dotted curves are the results for the P- and D-wave interference term, and the solid curves are the total results. 
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FIG. 12: The single-spin asymmetry in DIS production of charged pions off proton and deuteron targets, as function 

of X. The results are obtained using the parametrization of Ref. [120] for fi(x) at — 2.5 GeV^, and the light-cone CQM 
predictions for hlix) at the hadronic scale from Ref. [55, 77]. The dashed and dotted curves are obtained separating the 
contributions to hi{x) from S and P waves, respectively. The solid curves show the total results, sum of the S- and P-wave 
contributions. 




FIG. 13: The contribution from different angular momentum components to the single-spin asymmetry ^^"f-^'*'') in DIS 
production of charged pions off proton and deuteron targets, as function of x. The results are obtained using the light-cone 
CQM predictions for h^^l^'^{x) at the hadronic scale and the parametrization of Ref. [120] for fi{x) at = 2.5 GeV^. 

The dashed and dotted curves are obtained separating the contributions to h^^'^^'^i^) from the interference of S- and P-wave 
components and from the interference of P- and D-wave components, respectively. The solid curves correspond to the total 
results. 




FIG. 14: The contribution from different angular momentum components to the single-spin asymmetry ylj)'^'^*^'' '^^^ in DIS 
production of charged pions off proton and deuteron targets, as function of x. The results are obtained using the light-cone 
CQM predictions for /ijy^''(a;) at the hadronic scale and the parametrization of Ref. [120] for fi{x) at = 2.5 GeV^. 
The dashed and dotted curves are obtained separating the contributions to h^^^^'^{x) from the interference of L~ = -fl and 
Lz = — 1 angular momentum components, and from the interference of = and Lz = +2 angular momentum components, 
respectively. Solid curves: total results. 
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